Two-Loop Bethe Logarithms for non-5 Levels 
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Two-loop Bethe logarithms are calculated for excited P and D states in hydrogenlike systems, and estimates 
are presented for all states with higher angular momenta. These results complete our knowledge of the P and D 
energy levels in hydrogen at the order of a 8 m c c 2 , where m c is the electron mass and c is the speed of light, and 
scale as Z 6 , where Z is the nuclear charge number. Our analytic and numerical calculations are consistent with 
the complete absence of logarithmic terms of order (a/71-) 2 (Za) 6 \n[(Za)~ 2 ] m e c 2 for D states and all states 
with higher angular momenta. For higher excited P and D states, a number of poles from lower-lying levels 
have to subtracted in the numerical evaluation. We find that, surprisingly, the corrections of the "squared decay- 
rate type" are the numerically dominant contributions in the order (a/iv) 2 (Zctf m c c 2 for states with large 
angular momenta, and provide an estimate of the entire Beo -coefficient for Rydberg states with high angular 
momentum quantum numbers. Our results reach the predictive limits of the quantum electrodynamic theory of 
the Lamb shift. 
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I. INTRODUCTION 

The two-loop Bethe logarithm stems from a nonrelativistic 
treatment of the full two-loop self-energy. It is the finite part 
which is left over when the two-loop problem is renormal- 
ized according to the original derivation of Bethe [1] for the 
one-loop self-energy and therefore represents a natural gener- 
alization of the most basic quantum electrodynamic calcula- 
tion to the two-loop level, yet the respective calculation is the 
only way to complete the analysis of hydrogenic energy levels 
in the order (a/ir) 2 (Za) 6 m c c 2 . In comparison to the one- 
loop problem, the two-loop nonrelativistic self-energy is much 
more complicated, and it involves three instead of one propa- 
gator denominators. Matrix elements cannot be expressed in 
closed analytic form. Numerical methods, inspired by lattice 
calculations, represent convenient tools for the evaluations. 

Final calculations were carried out on workstation clusters 
of the Max-Planck-Institute for Nuclear Physics in Heidel- 
berg and of the National Institute of Standards and Technol- 
ogy. The total CPU time of all computations required for the 
present paper was about 30 months, but the real time was of 
course drastically shortened due to parallel processing. The 
two-loop problem is known to be computationally demanding. 
Additional difficulties arise for excited states due to a number 
of poles from lower-lying levels which have to be subtracted 
in the numerical evaluation. After subtraction, the result is fi- 
nite, as it should be, and represents an observable energy shift. 

The problem of so-called squared decay rates is associated 
with the two-loop Bethe logarithm [2]. The squared decay 
rates, which originate from specific photon energies where 
two propagators become singular simultaneously, cannot be 
uniquely interpreted as energy corrections even if the concept 
of a pole on the second sheet of the Riemann surface defin- 
ing the propagator is used in order to define energy levels, 
and have been assigned [2] to an contribution SBqq, where 
B$o is the two-loop coefficient multiplying the scaling fac- 
tor of (a/ir) 2 (Z a) 6 rUeC 2 / n 3 (here, n is the principal quan- 
tum number). For 3P and AP states, the squared decay rates 
were originally evaluated in in Ref. [2]. In [3], an analysis of 



S states was supplemented, and it has been clarified that the 
contributions as given in [2] should be understood in radians 
per second rather than cycles per second, or Hz. Here, we sup- 
ply numerical values for the correction of the "squared-decay" 
type for excited S, P, D, F, G and H states, and we also dis- 
cuss approximations for general states. The squared decay 
rates are a manifestation of the fact that a proper definition of 
energy levels ceases to be possible at the level of a 2 (Za) 6 in 
units of the electron mass, as already pointed out in Ref. [4]. 
Yet, these contributions are mathematically well defined and 
have to be evaluated in a calculation whose aim is to explore 
the predictive limits of the quantum electrodynamic theory of 
hydrogenlike systems. 

This paper is organized as follows. In Sec. II, we summa- 
rize all known results for the two-loop energy shift of excited 
P and D states up to the order a (Za) e (all energy shifts are 
measured in units of m c c 2 in this paper unless stated other- 
wise). Important definitions and formulas regarding the two- 
loop Bethe logarithm are recalled in Sec. EH. The asymptotics 
of the two-loop integrand and the calculation of the two-loop 
Bethe logarithm are discussed in Sec. IV. The squared decay 
rates, which are numerically significant especially for states 
with higher angular momenta, are treated in Sec. V. Miscel- 
laneous two-loop results are compiled in Appendixes A — C. 
One-photon vacuum-polarization and self-energy shifts of or- 
der a(Za) 7 are discussed in Appendix D. 



II. KNOWN TWO-LOOP RESULTS 

We work in natural units (h = c = eq = 1), as is customary 
in QED bound-state calculations, and we also set the electron 
mass equal to unity in the following. The real part of the en- 
ergy displacement of a hydrogenic state due to the two-loop 
(2L) self-energy (SE) can be parameterized as follows, 
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where H is a dimensionless function. Here, we are concerned 
with P and D states, and in part with states of higher angular 
momenta. For these manifolds, the first nonvanishing terms 
in the semi-analytic expansion of the dimensionless function 
H(Za) in powers for Za and ln(Za) read as follows, 



where k = (— lp'+'+V 2 (j + ±) is the Dirac angular quan- 
tum number (I is the orbital angular momentum quantum 
number, and j denotes the total angular momentum of the 
electron). In this way, one immediately reproduces the well- 
known general result 
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The first index of the ^-coefficients denotes the power of Za 
including the factors Za contained in Eq. (1), and the second 
index denotes the power of ln[(Za)~ 2 ]. For the lowest-order 
B40 term, a complete result is known which is valid for any 
hydrogenic state (see, e.g., Ref. [5]). We here express the re- 
sult in terms of the expectation value of an operator, as in 
Eq. (8.1) of Ref. [6], and suppress all operators which are zero 
for P states and states with higher angular momenta, 
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We adopt here all conventions of Ref. [6], in particular cr y = 
[a % , a^]/(2\) for the (2 x 2) -dimensional spin matrices. We 
recall the well-known identity 
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which is valid for all non-S 1 hydrogenic states. 

Our goal here is to extend formula (4) to two relative orders 
of Za, for general P and D states, and to find estimates for 
general states with higher angular momenta. In Ref. [6], nu- 
merical results have only been indicated for the fine-structure 
difference of D states [see Eqs. (4.22) and (5.3) of Ref. [6]] 
and for P states [see Eqs. (4.24) and (5.4), (6.2) and (7.7) of 
Ref. [6]]. However, the complete result for B 60 entails the 
two-loop Bethe logarithm, which has not been known so far. 

We recall that in terms of matrix elements, to be evaluated 
on Pauli-Schrodinger nonrelativistic wave functions of non-S 1 
hydrogenic states, the higher-order terms B^, Bqi and B^q 
can be expressed as follows (see Ref. [6]), 
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We take the opportunity to point out that in the corresponding 
Eq. (8.1) of Ref. [6], aprefactorof (a/7r) 2 in front of the sec- 
ond term on the right-hand side was missing [in Ref. [6], the 
entire above result (5) for the energy shift was multiplied by 
a factor of (a/ir) 2 on the left and right-hand side]. Also, the 
square of the logarithm In 2 [i(Za)~ 2 ] in the second and in 
the last-but-one term on the right-hand side were not included. 
Both typographical errors are absent from Eq. (8) of previous 
work reported in Ref. [7] and from the above Eq. (5). 

We recall here the definitions of the quantities N, (i^ and 
/?5, which enter into Eq. (5). In doing so, we first recall that 
in accordance with the notation introduced in [7], we redefine 
the finite part of an integral which diverges as 

A 

J dkf(k)=A\ + B\n\ + C + 0(\- 1 ), (6a) 
o 

with A = A(Za)~ 2 , for a specified upper cutoff A, to be 
equal to just the constant term, i.e. 



dkf(k)=C. 
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With this definition in mind, we have 
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The N term has previously been defined in Refs. [8, 9]; it is 
generated by a Dirac delta correction to the Bethe logarithm. 
The notation Sy is used in accordance with Ref. [6, 8] in order 
to denote the first-order perturbation of the Hamiltonian, the 
energy and the wave function in the ensuing matrix element, 
due to the specified potential V. 

The evaluation of the two-loop Bethe logarithm for 15 and 
25 has been discussed in Ref. [9], and for 35-65 in Ref. [3]. 
For 15 and 25, there is no ambiguity in the definition of the 
Bethe logarithm, because the integration over both photon en- 
ergies in the nonrelativistic self-energy is free of singulari- 
ties. However, for all higher excited 5 states and all states 
considered here, one incurs real (rather than imaginary) con- 
tributions to the energy shift from the product of imaginary 
contributions due to singularities along both photon integra- 
tions, and these result in "squared decay rates" in the sense of 
Ref. [2]. Thus, it is helpful to make a clear distinction between 
the singularity-free, principal-value part b^ and a real part 
S 2 Beo, which is incurred by "squared" (or, more precisely, 
products of) imaginary contributions from the pole terms. We 
write 



b L = b L + S 2 B t 
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where 6l is obtained as the nonlogarithmic energy shift stem- 
ming from the nonrelativistic self-energy, with all integrations 



carried out by principal value, and 5 2 B$o is the corresponding 
contribution defined in Refs. [2, 3], due to squared imaginary 
parts. The exact meaning of the separation (8) will be clari- 
fied below. Here, we just note that for 35-65 states, the above 
separation is not really essential, because S 2 Bqo is a numeri- 
cally marginal contribution as compared to \>l (see Ref. [3]), 
and thus b L {nS) « b L {nS) to a very good approximation. 
For the states under investigation here, the distinction (8), sur- 
prisingly, proves to be highly essential; yet before we come to 
a discussion of this surprising phenomenon, let us first discuss 
the evaluation of ~b~L- 



III. TWO-LOOP BETHE LOGARITHM 

We briefly recall [3, 9] the nonrelativistic two-photon self- 
energy A_Enrqed as an integral of the following structure, 

/ 2a \ 2 f £l f e2 
A^nrqed = f J / dkiki dk 2 k 2 /(fci, k 2 ) , 

(9) 

where the k\ and k 2 represent photon energies, and ei as well 
as e 2 are cutoff parameters. The function f(ki,k 2 ) is defined 
in Eq. (11) below; its precise structure is unimportant for the 
following consideration, which concerns the relation of the 
cutoff parameters e, (i — 1,2) to the ultraviolet cutoff A pa- 
rameter used in Eqs. (6) — (7). 

In order to clarify this relation, we recall that in the con- 
text of the e method (see §123 of Ref. [10] and [11-13]), 
the cutoff parameters are chosen so that the e can be made 
arbitrarily small, but only under the condition Za <C e, so 
that the expansion, first carried out in Za, then in e, for the 
low- as well as the high-energy parts, gives the complete re- 
sult for the self-energy. This procedure has been fully clari- 
fied in §123 of Ref. [10] and in Refs. [11-14]. In Ref. [14], 
it has been stressed that this method actually corresponds to 
an expansion in large e, and indeed, in the context of the di- 
mensional regularization method [6], the nonlogarithmic term 
which remains after subtraction of the divergent contributions 
as ei = Ai — ► oo, e 2 = A 2 — > oo, has been been identified 
as the two-loop Bethe logarithm 6l- Because we are dealing 
here with excited states that can decay via dipole radiation, 
care must be taken in the definition of the integration prescrip- 
tion for the k\ and k 2 integrations. In the current section, we 
first assume a principal-value prescription and write 
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where (P.V.) denotes the principal value, and it is understood 
that divergent terms for large Ai and A 2 have to be subtracted, 
in the sense of Eq. (6). 

The subtractions for large Aj (i = 1, 2) are carried out as- 
suming Ai <C A 2 , so we first let A 2 — ► oo, extract the con- 
stant term as a function of k\ and then we integrate this term 
with respect to k\, letting Ai — > oo. In this way we "sweep" 
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the entire first quadrant of the two-dimensional (fci , fo) plane 
(i.e., the entire region fci > 0, k 2 > 0). Of course, the same 
result would be obtained for the constant term under the op- 
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posite sequence of first letting Ai approach infinity, then A2. 

We are now in the position to recall the explicit form of 
f(ki, k 2 ), which reads 
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- ( 1 

+ 2 \ P H-E + k 2 



1 



1 



fci + k 2 \ H — E + k 2 



1 



1 



h + k 2 \ H - E + ki 



(11) 



r 



The interpretation of the terms on the right-hand side is as 
follows: the first six are due to fourth-order perturbation the- 
ory generated by the "velocity-gauge" nonrelativistic (p ■ A)- 
interaction. Herein, the fifth and the sixth terms are derivative 
terms which naturally occur in fourth-order perturbation the- 
ory. The seventh term involves (on the left and on the right) 
(p ■ ^-interactions, and a seagull term proportional to A 2 in 
between. The eighth and ninth terms involve a seagull term 
outside of the (p ■ A)-interactions. We also recall that the scal- 
ing of f(k\,k 2 ) with Za is completely clarified in Ref. [3]. 



IV. ASYMPTOTICS OF THE INTEGRAND 

We first consider the asymptotics relevant for the initial k 2 
integration in Eq. (10). In the limit k 2 3> fci, we find 
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The second term vanishes for D states and all states with 
higher angular momenta. Subtracting the two above terms ac- 
cording to the definition (6) [a practical procedure is outlined 
in Eq. (26) of Ref. [3]], we obtain 



5 (fci) = (P.V.) / dk 2 k 2 f(k u k 2 ). (13) 
Jo 



We now have to calculate 
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and this necessitates the calculation of the asymptotics of 
g(k\) for large fci. After a rather long and tedious calcula- 
tion, we find 
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a result which is valid only for non-S* states. The N term, 
which is defined in Eq. (7), remains important for the P state 
calculation, but vanishes for D states and all states of higher 
angular momenta. 

In the actual calculations, we use a lattice formulation of 
the Schrodinger propagator [15] with up to 200 000 lattice 
sites. This provides sufficient accuracy (and computational 
efficiency) for all calculations reported here and is feasible on 
average-size workstation clusters. 

It is computationally advantageous to isolate the contribu- 
tion to 6l due to the fourth term on the right-hand side of 
Eq. (8), in order to avoid numerical problems associated with 
the calculation of the reduced Green function. This contri- 
bution reads, after the fci and fc 2 -integrations and appropriate 
subtractions, 
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where <j> is the reference state. The calculation of this contri- 
bution is done separately, by writing it in terms of a sum over 
virtual states, which are calculated as a basis set on a lattice in 
coordinate space. Numerical values are compiled in Table I, 
where the name "loop-after-loop" should be taken cum grano 
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salis because the negative-energy virtual states are excluded 
from expression (16). The final results for b^ are in Table II. 

Based on the trend of the data in Table II, we conjecture 
here that in the limit n — > oo, the two-loop Bethe logarithms 
for a given I should approach a constant in the same way as the 
one-loop Bethe logarithms [16-18]. Also, the data in Table II 
indicate that the two-loop Bethe logarithms, just like their 
one-loop counterparts, become smaller in magnitude for in- 
creasing orbital angular momentum. Indeed, they do so quite 
drastically, with the D-state values being almost two orders 
of magnitude smaller than the P-state logarithms. The mag- 
nitude of the S-state logarithms [3], in turn, is in the range of 
60 ... 80 and thus larger than the P-state values by more than 
one order of magnitude. 



V. SQUARED DECAY RATES 



As clarified in Ref. [2] and Sec. IV of Ref. [3], squared de- 
cay rates represent natural limits to which energy levels can 
be uniquely associated with a particular atomic level. Cor- 
rections of this type have been analyzed for 2P and 3P in 
Ref. [2], and for 35 and AS in Ref. [3]. Here, we supplement 
values for 5P and 6P, as well as all other states with n < 6. 
Some inaccuracies in previous treatments for S and P states 
are corrected in Table III, and a more extensive list of levels is 
covered. 



TABLE I: Loop-after-loop contribution L 2 to the two-loop Bethe logarithm for P, D and F states. For the 
6H state, the result for the loop-after-loop contribution reads —0.000(1), implying a negative sign. 
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According to Refs. [2, 3], the term S 2 Bqq is generated by 
"squared decay rates," which correspond to well-defined, iso- 
lated points of the (fci, fc 2 )-photon energy plane, where two 
propagators become singular simultaneously in the integrand 
of the two-loop Bethe logarithm given in Eq. (11). These are 
all points where, for a given principal quantum number n of 
the reference state, any two of the following conditions are 
satisfied simultaneously, 
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where mi,m2,m < n are values of the principal quan- 
tum numbers in the intermediate states. In Fig. 1, we list 
all of these points in the two-dimensional (fci, fc2)-plane for 
a state with a principal quantum number n = 6. All fre- 
quencies which give rise to the singularities are smaller than 
(Za) 2 /2 = Z 2 R oc . 

Upon picking up the squared imaginary contribution, we 
obtain for the following structures from the two-loop inte- 
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One might ask why the result on the right-hand side reads 
-7r 2 , not -(27r) 2 . The answer is that both the fci and the fc 2 
integration contours have to be deformed along half circles, 
centered at the location of the singularities, infinitesimally be- 
low the real axis, i.e. in the mathematically positive sense. 
This half circle entails a factor 7ri, whose square gives the re- 
sults on the right-hand sides of Eqs. (20a) — (20d). 
Finally, the total squared decay rate 
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is obtained after picking up all terms which result as the prod- 
uct of two imaginary parts, when integrating over both fci and 
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FIG. 1: Location of the discrete pairs of frequencies (fci, £2) which 
give rise to the squared decay rates for an n = 6 state. The point 
Z 2 R oc — (Za) 2 m/2 denotes a "^-dependent Rydberg constant" 
and marks the ionization energy of a nonrelativistic hydrogenlike 
atom system. 



&2 in Eq. (11), according to the prescription (20). The con- 
tributions (5 2 i?60 is thus mathematically well-defined and, in 
particular, has an unambiguously defined sign for each indi- 
vidual state. 

A very curious observation can be made regarding the 
order-of-magnitude of versus that of S 2 Bqo for different 
manifolds of hydrogenic states. For S and P states, we have 
|(5 2 £?6o| <C |6i|, as is evident from Table II of this work and 
from Table II of Ref. [3]. For the IS and the 2S states consid- 
ered in Ref. [9], we even have 5 2 B eo (lS) = 5 2 B m {2S) = 
0. For higher excited S states, the magnitude of S 2 Bqo is 
smaller than the magnitude of the current numerical uncer- 
tainty of T>l incurred by the numerical integration described in 
Ref. [3], and thus, 5 2 Bqq can at present be completely ignored 
for S states. Therefore, the question of whether to include 
5 2 B 60 into the canonical definition of had been ignored in 
Refs. [3, 9] for the simple reasons that 5 2 B^(nS) was either 
vanishing or found to be a numerically negligible. Initially 
(see Ref. [2]), the notation S 2 B 60 was chosen such as to make 
an association with a small "uncertainty" in the definition of 
&£. However, the situation changes drastically for P states, 
where \S 2 Bqq\ ps 0.1 for n > 3, and even more drasti- 
cally for D states, where <5 2 -B60 3> for all n, revers- 
ing the hierarchy of the "uncertainty" <5 2 £> 6 o and the two-loop 
Bethe logarithm b^ (the latter would otherwise be assumed to 
constitute the dominant effect). 

In view of this situation, the pressing question arises 
whether or not one should include 5 2 B 6n (nS) into the defi- 



nition of the two-loop Bethe logarithm for states with higher 
angular momenta. One might argue that 5 2 B^ (nS) should be 
excluded from the definition of the two-loop Bethe logarithm 
bL because, as shown in Ref. [2], the term cannot be uniquely 
interpreted as an energy shift associated with a specific atomic 
level. On the other hand, one might argue that 5 2 B 60 (nS) 
should be included into the definition of the two-loop Bethe 
logarithm b^ for higher excited states, because it is generated 
by the product of two imaginary energy shifts, both of which 
originate from the same nonrelativistic two-loop self-energy 
(9), which otherwise gives rise to the low-energy contribution 
to the Bqo coefficient. We reemphasize that the <5 2 £>6o-term 
has an unambiguously defined sign for each individual state 
and is important for the comparison of our analytic approach 
to B eo to any nonperturbative (in Za) numerical calculation 
of the two-loop energy shift for the listed excited hydrogenic 
states, as well as for the comparison of theory and experiment 
based on a line-shape formalism [4]. Nonperturbative calcula- 
tions of the two-loop self-energy have recently been pursued 
intensively [19-25]. 

In the current investigation, we would like to follow the sec- 
ond route and define [see Eq. (8)] &l = &l + 5 2 Bq , where 
~b~L is obtained by performing all integrations in Eq. (10) via a 
principal-value prescription, with the results listed in Table II, 
and S 2 B & o contains all the terms generated by the products of 
imaginary parts incurred at poles in the (fci, fc2)-plane. This 
definition has the following advantages: (i) When the defini- 
tion (8) is adopted, the term b^ in the result for the B eo coef- 
ficient as given in Eq. (8. 1) of Ref. [6] contains all real (rather 
than imaginary) terms which can be inferred from the nonrel- 
ativistic self-energy given in Eq. (16) of Ref. [8] and Eq. (5) 
of Ref. [9] upon contour integration over both photon ener- 
gies, (ii) In nonperturbative numerical evaluations of radia- 
tive two-loop energy shifts for excited states, one would "pick 
up" terms of the type <5 2 i?60 naturally if one deforms both 
photon energy integration contours according to the Feynman 
prescription near adjacent bound-state poles. The definition 
(8) makes possible a direct comparison of nonperturbative 
numerical results obtained in this way, to the results for b^ 
and £>6o given here, without the need for any subtraction of 
the squared imaginary parts. (Hi) For D states, the definition 
(8) includes the numerically most significant real (rather than 
imaginary) contribution which can be inferred from the non- 
relativistic two-loop self-energy (9) into b^ and thus into the 
Bgo coefficient. 

Despite these advantages of the proposed definition (8), one 
should remember that the term 5 2 B eo cannot be interpreted as 
an energy shift uniquely associated with a particular atomic 
level, and has to be treated differently, namely according to 
ideas outlined in Sec. 3 of Ref. [2] when in a comparison to 
experiments. According to the discussion in the cited refer- 
ence, a valid way of treating (part of) 5 2 B^ would first entail 
a subtraction of this term from from the energy shift, and then, 
a reinterpretation of it in terms of an off-diagonal "decay rate 
operator" which has to be included into a line-shape formal- 
ism. 
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TABLE II: Numerical values for the principal- value integrated contribution 6l to the two-loop Bethe logarithms for P, D and F states. 
Values for S states are taken from [3] and are listed for completeness. For the 6H state, the result for b L reads —0.000(1), implying a 
negative sign. 



level 


b L 


level 


b L 


level 


b L 


level 


b L 


level 


b L 


level 


b L 


IS 


-81.4(3) 






















2S 


-66.6(3) 


2P 


-2.2(3) 


















3S 


-63.5(6) 


3P 


-2.5(3) 


3D 


-0.006(2) 














4S 


-61.8(8) 


4P 


-2.8(3) 


4D 


-0.004(2) 


4P 


-0.002(1) 










5S 


-60.6(8) 


5P 


-2.8(3) 


5D 


-0.005(3) 


5P 


-0.002(1) 


5G 


-0.001(1) 






6S 


-59.8(8) 


6P 


-2.9(3) 


6D 


-0.006(4) 


5P 


-0.002(1) 


6G 


-0.001(1) 


6H 


-0.000(1) 



TABLE III: Numerical values for the correction S 2 Bso for S, P, D and F states. Values for 35* and 4S have already been 
discussed in Ref. [3]. As explained in the text, the contributions listed below cannot be unambiguously interpreted as energy 
shifts associated with a particular atomic level, although they are mathematically well defined. 

level 8 2 B 6 o level <5 2 B 6 o level 8 2 B 60 level S 2 B 60 level <5 2 B 6 o level S 2 B 60 

~TS 00 =~ 

2S 0.0 2P 

3S -0.071 3P 

45 -0.109 4P 

55 -0.129 5P 

6S -0.141 6P 



-0.008 

-0.177 31) 0.130 

-0.243 4£> 0.183 4P 0.027 

-0.276 5D 0.203 5P 0.036 5G 0.009 

-0.295 6D 0.215 6P 0.044 6G 0.014 6H 0.004 



VI. EVALUATION OF B 60 



Before we give the general result for B 60 of P states, we would like to recall the known results for the other coefficients 
in Eq. (2). Specifically, we have for the spin-independent double-logarithmic term, B 62 (nP) = ■§f n ^ 1 (see Refs. [26, 27]). 
Furthermore, the logarithmic terms [6] are B 61 (nP 1/2 ) = § N(nP) + (±§ - ^ In 2) and B 61 (nP 3/2 ) = § N(nP) + 

In 2). Numerical values for N(nP) can be found in Eq. (17) of Ref. [28]. Based on Eq. (8.1)ofRef. [6] and 
on the results for obtained in this paper as well as standard analytic techniques for the evaluation of matrix elements, we are 
now in the position to give complete results for the B 60 coefficients of P states (see also Table IV), 



Beo(nP 1 / 2 ) = 



27517 209 1223 4 n 2 - 1 



■ + 



+ 



25920 288 n 960 n 2 27 n 



ln 2 (2) - 



38 n 2 - 1 
81 n 2 



ln(2) 



25 

IT 



3 

2n 



9 
2^2 



C(2)ln(2) 



9151 



1 



+ ■ 



1009 



10800 An 1200 n 2 



C(2) 



25 
"24 



3 

8n 



9 
8^2 



C(3) 



B e o(nP 3 / 2 ) 



+ 04(nP 1/2 ) + /? 5 (nP 1/2 ) + 
73321 185 



38 4 , ^ 

ln(2 

45 3 W 



8111 



N(nP) + b L (nP) 
4 n 2 - 1 



ln 2 (2) 



11 n 2 - 1 
81 n 2 



ln(2) 



(22a) 



+ 



53 



/299 3 
V"80" ~ 8^ ~ 20n 2 



C(2)In(2) 



24377 1 

+ 



3187 



_299 _3_ 53 
21600 ' 16 n 3600n 2 I ^ ' + \ 320 + 32n + 80n 2 



C(3) 



■ MnP 3 /2) + MnP 3/2 ) + 



38 4 , . . 

ln(2 

45 3 W 



N(nP) + b L (nP) . 



(22b) 



We see that the total value of Bqq is the sum of high-energy operators given by inverse powers of the principal quantum number, 
linear combinations of terms proportional to ln(2) and £ functions of various arguments, and low-energy terms /J4 and /3 5 which 
are known from one-loop calculations [28, 29] (see also Tables VII and VIII below), as well as the two-loop Bethe logarithm b L . 

For D states, it is known that Bq 2 (jiD) = BQi(nD) = 0. Indeed, the fact that N(nD) — for D states implies that the 
term of order 1 jk\ in the asymptotics (15) vanishes, and this is consistent with the zero result for B &1 (nD). The nonlogarithmic 
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terms read [note in particular N(nD) = 0], 



125863 1021 18811 /1387 9 313 > 
Beo(nD 3/2 ) = -^ooo " 9600^ + 60480^ + 1,8400 + 40^ " 420^ ) C(2) ln(2) (23a) 



( 5501 3 1073 \ / 1387 9 313 \ A ,„, „ . „ . „ , „ , ? . . 

+ (" 151200 - 80^ + 7560^ j C(2) + {' 33600 " 160^ + 1680^ j C(3) + + + ^ ' 

„ . „ . 61133 949 4967 / 43 1 61 \ , ,„„, , 

^o(^5/2) = + - 30240^ +1,3150-10^+210^; C(2)ln(2) (23b) 

+ (" sloO + (L - 94^) « 2 > + (-12^0 + 40^ - 84^) C(3) + + A ^ + ^ ' 



r 



TABLE IV: Values of the total Beo coefficient for P states, derived 
from Eq. (22). For the numerical values of f3i and /3s, see Tables VIII 
and IX in Appendix A. 



level 


^60 


level 


#60 


2-Pl/2 


-1.6(3) 


2P3/2 


-1.8(3) 


3Pl/2 


-2.0(3) 


3P 3 /2 


-2.2(3) 


4Pi/ 2 


-2.4(3) 


4P 3 /2 


-2.5(3) 


5Pl/2 


-2.4(3) 


5P 3/2 


-2.5(3) 


6P1/2 


-2.5(3) 


6P3/2 


-2.6(3) 



TABLE V: Same as Table IV, but for D states. The numerical val- 
ues of /?4 and /3s, which are needed for the numerical evaluation 
of the expressions in Eq. (23), are recorded in Tables IX and X in 



level 


P60 


level 


^60 


3O3/2 


0.141(2) 


30 5/2 


0.123(2) 


4£> 3/2 


0.199(2) 


4A./2 


0.178(2) 


5£> 3/2 


0.219(3) 


5A/2 


0.196(3) 


6-D3/2 


0.230(4) 


6D5/2 


0.207(4) 



For F states, the values of B 6a are displayed in Table VI (de- 
tailed formulas are given in Appendix C). We observe that the 
total value of B ea for F states is numerically rather close to 
the value of <5 2 £>60 for each state as given in Table II. Con- 
tributions from the fine-structure dependent terms as well as 
those from the high-energy operators gradually vanish as the 
angular momentum increases, and the dominant remaining 
contribution then stems from S 2 Beo- 

For G and H states, the total value of B 60 is the same as 
the value of 5 2 B 60 up to the level of ±1.0 x 10 ~ 3 in units of 
B 60 , as shown in Tables III and VI. We may now use our ex- 
perience regarding the asymptotic behavior of one-loop Bethe 
logarithms [18] and of relativistic Bethe logarithms [29], and 
based on the data in Tables IV, V and Table VI for an extrapo- 
lation of £> 60 to an arbitrary state. Specifically, we conjecture 
that for a given hydrogenic energy level which reads nLj in 
the usual spectroscopic notation, 

B 6 o(nLj) « S 2 B 60 « (l.7 - ^ _ 2 £ ° + J ± 50% , 

for L > 3, and both j = L ± 1/2 . (24) 

This conjecture implies that the dominant contribution to B^q 
comes from the squared decay rates, i.e. from the term S 2 B eo , 
and the conjecture (24) can be used to estimate this coefficient 
for an arbitrary Rydberg state of high angular momentum. 



TABLE VI: Same as Tables IV and V, but for F, G and H states. 
For these states, the total values of B60 are close to the numerical 
values for 8 2 Be,o as listed in Table III, and the high-energy operators 
as well as the fine-structure dependent corrections do not introduce 
any significant numerical deviation of the total value of B60 from 
S 2 B 60 - 



level 


#60 


level 


P60 


4P 5 /2 
5P 5 /2 
6F5/2 


0.027(1) 
0.037(1) 
0.046(1) 


4P 7 /2 
5P 7 /2 
6P7/2 


0.023(1) 
0.033(1) 
0.041(1) 


5G7/2 
6G 7/2 


0.009(1) 
0.014(1) 


5G9/2 
6G9/2 


0.008(1) 
0.013(1) 


6-Hg/ 2 


0.004(1) 


6-Hri/2 


0.004(1) 



VII. CONCLUSIONS 

Together with the one-loop results from Appendix D, the 
reported calculation of the two-loop nonlogarithmic term of 
order (a/ir) 2 (Za) 2 m c c 2 completes the analysis of the quan- 
tum electrodynamic corrections to non-S* states in the order 
a 8 m c c 2 for hydrogen (Z = 1). Together with other recent 
theoretical investigations [6, 7], the current work clarifies fur- 
ther prerequisites and the theoretical basis for the determina- 
tion of fundamental constants from hydrogen spectroscopy at 
the level of one part in 10 14 . For ionized helium, the correc- 
tions calculated here are enhanced by a factor of 2 6 = 64 in 
frequency units. The numerical values of principal-value two- 
loop Bethe logarithms b L are given in Table II, and squared 
decay rates S 2 Bqq are indicated in Table III. Final numeri- 
cal results for B eo , which is the sum of ~El, S 2 B &0 and high- 
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energy operators as well as fine-structure dependent terms 
[see Eqs. (22a) — (23b)], are given in Tables IV, V and VI. Our 
total results for Bqo contain the sum of the contributions from 
all four gauge-invariant subsets as shown in Fig. 2. We con- 
firm that the trend already observed in the literature [17, 18] 
for one-loop Bethe logarithms which decrease in magnitude 
with increasing orbital angular momentum quantum number, 
also holds for the two-loop counterparts. Finally, a com- 
parison of our results with a very recent numerical inves- 
tigation [25] indicates that our results for i? 6 o(2P 1 / 2 ) and 
S 6 o(2P 3 / 2 ) for subset i (see Fig. 2 and Tables VII and VIII 
below) are consistent with a nonperturbative (in Za) treat- 
ment of the two-loop self-energy, which forms subset i in the 
classification of Fig. 2. 



o 







our calculation explores the predictive limits of quantum elec- 
trodynamic theory of energy levels of the hydrogen atom and 
hydrogenlike systems. Explicit theoretical values of the Lamb 
shift for selected states in hydrogenlike ions will be presented 
elsewhere. 

All calculations reported here (see in particular Tables II 
and III) were carried out on an average-size cluster of work- 
stations. A fifth-order finite-difference scheme was used in or- 
der to accurately describe the behavior of the S component of 
the propagators on a grid near the origin, and in order to allow 
for a fast calculation of the two-loop integrand on the work- 
station cluster of Max-Planck-Institute for Nuclear Physics 
in Heidelberg and of the National Institute of Standards and 
Technology. However, even the most sophisticated computer 
hardware and numerical algorithms would be useless in the 
current context unless combined with a thorough formulation 
of the subtraction procedures necessary to extract the physi- 
cally observable energy corrections and frequency shifts. 



o 

Q 



FIG. 2: The two-loop diagrams which give rise to the Beo coefficient 
naturally fall into four gauge-invariant subsets, labeled i-iv (we fol- 
low the classification of diagrams according to Ref. [6]). Our results 
for B 60 contain contributions from all of these diagrams; the sepa- 
ration into the four subsets is only relevant for the detailed formulas 
given in the Appendixes. 



The connection of the two-loop problem to squared decay 
rates has been analyzed previously in Refs. [2, 3]. For the 
atomic levels under investigation, we give a respective treat- 
ment in Sec. V. The contributions to B 60 in Table III due 
to squared imaginary parts cannot be interpreted as energy 
shifts in the usual sense and set a limit to the actual "defin- 
ability" of the atomic energy levels under investigation. The 
surprising conclusion of the current investigation is that the 
squared decay rates actually give the dominant contribution to 
£?6o for states with higher angular momenta (see Tables II, III 
and VI). As S 2 B &0 cannot be associated uniquely as an en- 
ergy correction to one specific atomic level, we conclude that 
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APPENDIX A: MISCELLANEOUS RESULTS FOR P 
STATES 

In Sec. VI, we had discussed the evaluation of the total Bqq 
coefficient for P 1 / 2 and P 3 / 2 states of general n. Here, we 
would like to supplement the results for the individual sub- 
sets i-iv of the two-loop diagrams as given in Figs. 1-4 of 
Ref. [6] and in Fig. 2 here. The results presented in this Ap- 
pendix might be helpful both for a comparison of our analytic 
to nonperturbative numerical results, as well as for a verifica- 
tion of alternative, independent analytic calculations for par- 
ticular sets of diagrams. 

We start with the subset i (pure two-loop self-energy) for 
Pi/2 states (see Fig. 2), 



Bl {nPy 2 ) 



3829 89 



44437 



4 n 2 



1 . 2 ,„, 38 n 2 - 1 



77760 96 n 77760 n 2 27 n 2 

7613 



ln'(2)- 



81 n 2 



ln(2) 



Hf^-^K(2)ln(2) + 



5 1663 \ / 25 

5400 4^ + 900n 2 ) ^ ' + \ 24 



+ f34nP 1/2 ) + f3 5 (nP 1/2 ) + 



10 



ln(2) 



N(nP) + b L (nP). 



(Al) 
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TABLE VII: Breakdown of the individual contributions to S6o(nPi/ 2 ) from all diagrammatic subsets as de- 
fined in Ref. [6]. The numerical uncertainty in Be () (nP 1 ^ 2 ) is entirely due to the two-loop Bethe logarithm 
given in Table II. The N term has been calculated in Ref. [28]; it is spin-independent and therefore indicated 
only for the P 1/2 states (cf. Table VIII below). 

level 0* p 5 N B^o Bffi B$ fljg 

2P 1/2 -0.314 100 0.260 017 0.003 301 -1.5(3) 0.017 642 -0.034 214 -0.126 543 

3Pi /2 -0.362 439 0.292 936 0.003 572 -1.8(3) 0.022 297 -0.040 459 -0.149 977 

4P 1/2 -0.380 983 0.304 244 -0.000 394 -2.2(3) 0.024 123 -0.041 562 -0.158 179 

5Pi /2 -0.390 543 0.309 438 -0.004 304 -2.2(3) 0.025 057 -0.041 519 -0.161 975 

6Pi/ 2 -0.396 287 0.312 249 -0.007 497 -2.3(3) 0.025 613 -0.041 211 -0.164 037 



We also note that in Eqs. (5.4), (6.2) and (7.7) of Ref. [6], results for the fine-structure difference associated with the diagrams 
of subsets ii-iv were indicated. We reemphasize that these results are valid only for the fine-structure difference of P states. 
Therefore, we would like to supplement the individual contributions for P 1 / 2 and P 3 / 2 here, starting with P1/2, 

*eo(nP 1/2 ) = - ^ - ^ + ^ + + \ - ^ C(2) , (A2) 



1939 


119 


2160 


72 n 






-{^ 




45 V 





B»P 1/2 ) = - |(l - -L) - ± N{ nP) , < Kl/2 ) = (l - -1) . (A3) 

For the two-loop self-energy subset i, we obtain in the case of a P 3 / 2 state, 

„, . „ , 73321 185 8111 4 n 2 - 1 , 2 ., 11 n 2 - 1 , 

Bm(nP 3/2 ) = 1 1 ^ H ^— n (2 ^— ln(2) 

60V 3/2^ 10368O 1152 n 25920 n 2 27 n 2 W 81 n 2 y ' 

/299 3 53 \ / 24377 1 3187 \ / 299 3 53 \ 

+ (^C7 - 8^ - 20^j C(2)M2) + (-21600 + 16^- 3600^J C(2) + {' 320 + 32^ + 80^ j C(3) 



+ l3i(nP 3/2 ) + /3 5 (nP 3/2 ) + 



N(nP) +b L (nP) . (A4) 



The contributions of the diagrammatic subsets ii-iv read as follows, 

3595 119 11 (Yi 
5184 + 288n ~ 88n 2 + I 40 An ' 90n 2 



^oKs /2 ) - -S + ^- ^+f^-i + F^) C(2) 



^(nP 3/2 ) = - 1 (l - ^) - 1 N(nP) , ^ (nP 3/2 ) = - ^ (l - ^) . (A6) 

It is perhaps also instructive to indicate here, the total result for the fine-structure difference of P 6 o- ln Ref- [6], only the 
individual subsets i-iv were treated for the fine-structure, but their sum was not indicated, and this "checksum" is supplemented 
here: 

/ 7-j \ „ , „ s 1361 1021 2491 ln 2 -l / 103 15 37 \ „ nl , 

B eo (nP 3/2 ) - B 6o{ nP 1/2 ) = — + — - — 2 + - — In(2) + (-^ - ^ + ^ j C(2) ln(2) (A7) 

/ q s 2 \ / 1 03 15 37 \ 

+ "32 + TTn + 45^ + 960 + 3^ " 80^ ^ + ~ + ~ " 



In this result, 6l cancels, and P6o( n ^3/2) — B^(nP 1 / 2 ) evaluates, e.g., to a numerical value of —0.134817 for n = 2. The 
reader may consult Tables VII and VIII for details and observe that the principal numerical uncertainty due to b L cancels for the 
fine-structure. 
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TABLE VIII: Same as Table VII, but for nP 3/2 . 



level 


At 




Bio 


#60 


Tiiii 
#60 


#60 



2P 3/2 0.157 050 -0.130 009 -1.8(3) 0.004 207 0.015 786 0.063 272 

3P 3/2 0.181 219 -0.146 468 -2.2(3) 0.005 208 0.018 801 0.074 989 

4P 3/2 0.190 492 -0.152 122 -2.6(3) 0.005 510 0.020 938 0.079 090 

5P 3/2 0.195 271 -0.154719 -2.6(3) 0.005 627 0.022 481 0.080 988 

6P 3/2 0.198 144 -0.156 124 -2.7(3) 0.005 678 0.023 604 0.082 019 



APPENDIX B: MISCELLANEOUS RESULTS FOR D STATES 



For D states, we observe that only subsets i and ii contribute and start with D 3 / 2 , 

i 405431 453 4901 /1387 9 313 \ 

B 60 (nD 3/2 ) - -^ooo + 320TM " 12096 n 2 + 1,8400 + 40^ " 420^ ) C(2) M2) 

/ 3469 3 4349 \ / 1387 9 313 \ 

+ V 30240 ~ 16n + 7560 n 2 ) ^ ' + \ 33600 ~ 160 n + 1680 n 2 ) **' ' 

+ 0i(nD 3/2 ) + (3 5 (nD 3/2 ) + b L (nD) . (Bl) 

Subset ii yields, 

„ u . „ , 5593 119 1547 / 47 3 13 \ 

BUnD m) = -— - — + ^ + ^ + ^ - ^ C(2) • (B2) 



We now continue with D 5 / 2 , 
21503 

756000 800 n ' 6048 n 2 ' V 3150 10 n ' 210 n 2 



i 21503 53 1577 / 43 1 61 

B 60 {nD 5/2 ) = -37777^7: - 77777— + 7777777-3 + tttttt, ~ 77^7 + 77777^2 C(2) ln(2) 



+ (Ho + IS - 5s£?) <(2 » + ("iio + 557. - sith) a3 > + M" D >r-> + A(»^») + «« p > ■ "»> 



„ w , \ 1819 119 17/107 1 9 , 

BUnD h/2 ) = —777 + 7777777- - 777-2 + - 77-7 + ^~ 2 ) C(2) ■ (B4) 



Here, subset m yields 

RiU™ n. ,„) = 

48600 1080 n 40 n 2 V 4725 15 " 35 
Subsets m and iv do not contribute in either case for D states. A "checksum" for the total fine-structure difference can be useful, 

„ , „ , „ , „ \ 777473 2597 5749 / 3817 13 29 \ , 
i? 60 (nZ? 5/2 ) - S 60 (^ 3/2 ) = TTjggg^gQ + 1^80^ " 1^096^ + \~ 25200 " 40^ + 28^ J C(2) M2) 

3817 13 29 \ / 3817 13 29 \ 

151200 + 240 n ~ 168 n 2 ) **' ' + ^100800 + 160 n ~ 112 n 2 ) ' 

+ MnD 5/2 ) - fa(nD 3/2 ) + f3 5 (nD 5/2 ) - (3 5 (nD 3/2 ) . (B5) 

In this result, cancels, and the quantity B 60 (nD 5 / 2 ) — B^{nD 3 / 2 ) evaluates, e.g., to a numerical value of —0.018197 for 
n = 2. The interested reader may consult Tables IX and X for details, observing that the principal numerical uncertainty due to 
bh cancels for the fine-structure. 
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TABLE IX: Breakdown of the individual contributions to Beo {nD 3 / 2 ) from all diagrammatic subsets as defined 
in Ref. [6]. The numerical uncertainty in Bl (nD 3 / 2 ) is entirely due to the two-loop Bethe logarithm given in 
Table II. 



level 


ft 


ft 


Bio 


&60 


3D 3 /2 

4£> 3/2 

5£>3/2 
6-D3/2 


-0.002 361 
-0.002 883 
-0.003 101 
-0.003 200 


0.005 397 
0.006 280 
0.006 675 
0.006 886 


0.141(2) 
0.199(2) 
0.220(3) 
0.230(4) 


-0.000 629 
-0.000 696 
-0.000 714 
-0.000 716 



TABLE X: Same as Table IX, but for nD 5/2 . 



level 


ft 


ft 


Bio 


r>ii 
-t>60 


3D 5 /2 
4£> 5/2 

5£> 5 /2 

6D5/2 


0.001 574 
0.001 922 
0.002 067 
0.002 133 


-0.003 598 
-0.004 186 
-0.004 450 
-0.004 591 


0.123(2) 
0.177(2) 
0.196(3) 
0.207(4) 


0.000 128 
0.000 182 
0.000 202 
0.000 209 



APPENDIX C: MISCELLANEOUS RESULTS FOR F, G AND H STATES 

For states with angular momenta I = 3, 4, 5, we indicate here only the final results for Bqo, without considering the breakdown 
for the terms generated by the individual subsets in Fig. 2. We obtain for F 5 j 2 states, 

8321 1415 2567 / 71 1 209 \ 
60 5/2) - - 666792 - 42336 - + lgl44 n 2 + ^2205 + Tin ~ 630 n 2 J ^ ^ n ^ 

/ 2173 1 347 \ / 71 1 209 \ 

+ V 317520 ~ 84n + 5670n 2 J ^ ' + \ 8820 ~ 56n + 2520n 2 J C( ' 

+ Pi(nF 5/2 ) + p 5 (nF 5/2 ) + b L (nF) . (CI) 

For F 7 / 2 states, the result reads 

„ , N 994501 1343 33031 / 125 9 89 \ „ „ , , 

B m (nF 7l2 ) = + - + + - C 2 In 2 

buv 7/2; 284497920 75264n 362880n 2 V 56448 224n 504n 2 / K ' 

5629 ' 3 1067 V(2) + f^ + ^ _ 89 A c(3) 



5080320 448n 45360n 2 /^ w V 225792 896n 2016n 2 / 
+ Pi(nF r/2 ) + fc(nF 7/2 ) + b L {nF) . (C2) 

G 7 j 2 states provide us with the following information, 

762011 67 2389 / 2677 1 1037 \ 

6o(n 7/2) - - Y91600640 ~ 4608 n + 29568 n 2 + V 266112 + 32^ ~ 5544 n 2 ) ^ > n ^ > 

16601 1 + ^W)+(-^-^ + ^W) 



7983360 192 n 99792 n 2 / ^ w V 1064448 128 n 22176 n 2 / 

ft(«G 7 / 2 ) + ft(nG 7/2 ) + b L (nG) , (C3) 
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whereas Gg / 2 states yield 



170587 193 1277 / 44 1 73 \ 

60 9/2) - 112266OOO + 21600 n ~ 22176 n 2 + V 23625 ~ 50^ + 630 n 2 J ^ ' n{ - > 

( 1153 1 1037 ^ r(9) ( 11 1 73 ^ rcn 

+ 12474000 + 300 n ~ 62370 n 2 ) ^ ' + \23625 + 200 n ~ 2520 n 2 ) ^ ' 

+ MnG 9/2 ) + /3 5 (nG 9/2 ) + b L (nG) . (C4) 

The results for H g j 2 

9169711 2203 19361 / 14411 9 1543 \ 

eo(n 9/2) - - 5606172 000 _ 290400 n + 370656 n 2 + \ 3539250 + 550 n ~ 12870 n 2 J ^ > n{ - > 

104891 3 1241 \ , , / 14411 9 1543 \ 

C(2)+ - 777^7^ -^773 + 777^2 C(3) 



127413000 1100 n 57915 n 2 J bW V 14157000 2200 n 51480 n 2 , 
+ fa(nH 9/2 ) + fc(nH 9/2 ) + b L (nH) , (C5) 

and fl"ii/ 2 

4047937 2131 29401 / 85 1 419 \ 

B 6Q (nH n/2 ) - 5381925120 + iT ^ T ^ - 741312n2 + ^75504 " 88^ + 5148^ ) C(2) ln(2) 

/ 877 1 1123 \ / 85 1 419 \ 

+ V 20386080 + 528 n ~ 92664 n 2 ) ^ > + ^302016 + 352 n ~ 20592 n 2 ) ^ ' 

+ (1 4 (nH u/2 ) + f3 5 (nH n/2 ) +b L (nH) (C6) 
complete the investigation of B 60 coefficients. 



APPENDIX D: VACUUM POLARIZATION AND SELF 
ENERGY SHIFTS OF ORDER a (Za) 7 

As is well known, the one-loop (1L) vacuum polarization 
(VP) correction scales as 

A^-^^vp^), (Dl) 

where Fvp(Za) is a dimensionless function. The first few 
terms of Fvp(Za) for S states can be found in a number of 
review articles, e.g. [30, 31]. E.g., the lowest-order term for 
S states is Aj^(nS) = -4/15. The first index of the A YP - 
coefficients denotes the power of Za, and the second denotes 
the power of ln[(Za:)~ 2 ]. 

For states with nonvanishing angular momenta, it is well 
known that the leading terms in the expansion of Fyp(Za) 
vanish. In the order a(Za) e , the leading coefficients for P 
states [32] are known to read as follows, 

A^{nP 1/2 )= (D2a) 
A^(nP 3/2 )=-^-^±. (D2b) 



An investigation of the behavior of the wave functions near 
the nucleus, inspired by Ref. [33], leads to the following cor- 
rection terms of order a(Za) 7 , 

These results are consistent with the particular case n = 2 
treated on p. 124 of Ref. [31], for which the coefficients read 
A 70 (2P 1/2 ) = 4171-/3072 and A 70 (2P 3/2 ) = 7tt/1024. Both 
the Aq P as well as the A^q coefficients vanish for D states 
and all states with higher angular momenta. Note that the self- 
energy remainder function Gse, as calculated in Refs. [34, 35] 
for P states, contains all contributions of order a {Za) 7 due 
to the one-photon self-energy. For D states and states with 
higher angular momenta, self-energy shifts of order a (Za) 6 
have been compiled, e.g., in Ref. [14], and the self-energy cor- 
rections of order a (Za) 7 vanish. These remarks supplement 
the above results for the one-loop vacuum polarization. 
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